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ABSTRACT 

The coalescence of a supermassive black hole binary (SMBHB) is thought to be 
accompanied by an electromagnetic (EM) afterglow, produced by the viscous infall of 
the surrounding circumbinary gas disk after the merger. It has been proposed that once 
the merger has been detected in gravitational waves (GWs) by the Laser Interferome- 
ter Space Antennae (LISA), follow-up EM observations can search for this afterglow 
and thus help identify the EM counterpart of the LISA source. Here we study whether 
the afterglows may be sufficiently bright and numerous to be detectable in EM surveys 
alone. The viscous afterglow is characterized by an initially rapid increase in both the 
bolometric luminosity and in the spectral hardness of the source. For binaries with 
a total mass of 10^ - IO^^Mq, this phase can last for years to decades, and if quasar 
activity is triggered by the same major galaxy mergers that produce SMBHBs, then it 
could be interpreted as the birth of a quasar. Using an idealized model for the post- 
merger viscous spreading of the circumbinary disk and the resulting light curve, and 
using the observed luminosity function of quasars as a proxy for the SMBHB merger 
rate, we delineate the survey requirements for identifying such birthing quasars. If 
circumbinary disks have a high disk surface density and viscosity, an all-sky soft X- 
ray survey with a sensitivity of Fx ^ 3 x 10"^"^ erg s"' cm"^ which maps the full sky 
at least once per several months, could identify a few dozen birthing quasars with a 
brightening rate dlnFx/dt > 10% yr"' maintained for at least several years. If > 1% 
of the X-ray emission is reprocessed into optical frequencies, several dozen birthing 
quasars could also be identified in optical transient surveys, such as the Large Syn- 
optic Survey Telescope. Distinguishing a birthing quasar from other variable sources 
may be facilitated by the monotonic hardening of its spectrum, but will likely remain 
challenging. This reinforces the notion that observational strategies based on joint 
EM-plus-GW measurements offer the best prospects for the successful identification 
of the EM signatures of SMBHB mergers. 

Subject headings: accretion, accretion disks — black hole physics — gravitational 
waves — quasars: general 



'Kavli Institute for Theoretical Physics, UCSB, Santa Barbara, CA 93106 



-2- 



1. Introduction 



Observational evidence robustly indicate s that all or nearly all galaxies harbor a supermas- 
sive black hole in their nucleus (SMBH; e.g.. ^Mago rrian et aL 1998 ). Since cosmological struc- 
ture formation models predict a hierarchy of galaxy mergers, if nuclear SMBHs were indeed 
common at earlier times, then th ese mergers should result in the formation of SMBH binaries 
(SM BHBs; feegelman et al]ll980l) and these binaries should then be common throughout cosmic 



time (Haehnelt 1994 :lMenou et al.ll200ll:IVolonteri et al.ll2003l : IWvithe & Loebll2003l : ISesana et al 
2007bl : lLippai et al.ll2009l : iTanaka & Raima J2009I) . 



It has also long been known , both observationally (e.g. JSanders et al.lll988h and theoretically 
(e.g., iBames & Hemquistlll99l|) that galaxy mergers can drive gas to the nucleus of the merger 
remnant, which could facilitate the merger of the nuclear SMBHs on one hand, while also provid- 
ing fuel for quasar activity on the other. Mergers are therefore generically also believed to trigger 
quasar activity; the rate of major galaxy mergers can indeed provide an explanation for the ob- 
serv ed evolution of the qu asa r population as a who le (ICarlberglll99d : for more recent work, see, 
e.g., iHopkins et al.ll2007al and lWyithe & Loebll2009l and references therein). 



Despite their expected ubiquity, observational evidence for SMBHBs is scarce, and the pre- 
cise timing of any quasar activity, and, when it occurs relative to the merger of the nuclear SMBH 
binary, remains unclear ([Kocsis et al.l |2006|) . A handful of pai rs of active SMBHs in the same 
galaxy have been resolv ed directly, at ~kpc separation in X-ray (iKomossa et al.l 2003) and o ptical 
(|Comerford et al. I I2OO9I) images, and at ~10pc separation in the radio (Rodriguez et al.ll2006|) . con- 
firming that gas is present around the SMBH binary, and that quasar activity can, at least in some 
systems, commence prior to their coalescence. However, there has been at least one suggestion 
that luminous activity can be occurring later, at the time of the merger, a s well - moment arily in- 
terrupted by the coalescence of the SMBHs and reactivated after-wards (iLiu et al.ll2003bl). While 



there are many more obse r ved SMBHB candidat e s with small separations (e.g 



Schoenmakers et al.ll2000l: iMerritt & Ekersll2002l : ISudou et all 120031: iLiul 12004 iBoroson & Lauer 



Roos et al. 1993 



20091), the evidence for these tighter binaries is indirect, and each candidate system has alterna- 



tive explanations. The expectation is that at large separations, the binaries rapidly lose orbital 
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compact, gravitational wave (GW) emission rapidly shrinks the orbit, culminating in a merger. 
How long this process lasts, and at what stage(s) the SMBHs light up as luminous quasars, is, 
however, also poorly understood theoretically. 



Apart from the cosmological context, interest in EM signatures of SMBH mergers surged re- 



cently (e.g., 
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driven by (i) the prospect that the Laser Interferometer Space Antenna^ {LISA) w ill detect the 



merg ers in GWs and provide a tractable list of (perhaps as low as a few hundred; e.g.. lKocsis et al. 



2008) EM candidates for SMB HBs and (ii) the breakthrough in numerical general relativi stic cal- 
culations of BH mergers (e.g.. |Pretoriusll2005l ; ICampanelli et al.ll2006l ; iBaker et al.ll2006|) . which 
led to robust predictions of significant mass-loss and recoil that can significantly perturb the am- 
bient gas. A simultaneous observation of the merger in gravitational and EM waves would enable 



new 



s cientific investigations in cosmology and BH accretion physics ( Cutlej 1998 



20051 ; iKocsis et al.ll2006L 120071 ; iLang & Hughesll2008l ; |Phinnevll2009l ; iBloom et al. 



-lolz & Hughes 



2009) 



In this paper, we focus on one particular signature of SMBHB coalescence, whic h we will 



hereafte r refer to as the "viscous afterglow". The physics of this model was discussed by lLiu et al. 



(l2003bl) in the context o f the interruption of jets in double-double radio galaxies, and later by 
Milosavljevic & Phinneyl (2005; hereafter MP05) in the context of EM counterparts of LISA sources. 
Prior to merg er, the SMBHB torques open and maintain a cavity in the center of a thin circumbi- 
nary gas disk (|Artymowicz & Lubowll 19941) . When the binary becomes sufficiently compact, GW 
emission causes the binary orbit to shrink faster than the gas just outside the cavity can viscously 
respond. The merger takes place inside the cavity, which is subsequently filled as the disk viscously 
spreads inward. Because the refilling inner disk produces higher-energy photons than the outer re- 
gions, the disk is predicted to transition from an X-ray-dim state to an X-ray-bright one, with its 
bolometric luminosity increasing by a factor of ~ 10 during this time. This transition is expected to 
take place on humanly tractable timescales, with the cavity filling in ~ 10(1 + z){M / \O^Mp,y-^ yr, 



where M is the total mass of the binary. A study of an optically selected sample by lOibson et al. 



{20081) found X-ray-dim AGN to be rare (^ 2% at z ~ 2), suggesting that it would be tractable to 
catalog and monitor such systems for possible observational signatures of a merger afterglow. 

In the observational scenario originally proposed by MP05, LISA would detect the GWs from 
the merger and determine its approximate location in the sky to within ~ 0. 1 deg, triggering a 
follow-up search to identify the EM counterpart and host galaxy. A natural question to ask, how- 
ever - and the subject of the present paper - is whether the viscous afterglows may be sufficiently 
bright and numerous to be detectable in EM surveys alone, even before LISA is launched. The 
identification of mergers by their EM signatures alone could, in fact, be valuable for several rea- 
sons. First, LISA will be sensitive to GWs from relatively low-mass SMBHBs, with total masses 
of ~ (lO'^- 10')/(1 -l-z)M0. EM studies could, in principle, detect coalescing SMBHBs outside 
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this mass range, and therefore complement the LISA binary population. Second, while many mod- 
els for the cosmological evolution of SMBHs predict that LISA wi ll detect dozens or h undreds of 
mergers (if "seed" black holes are abundant and merge often; e.g.. ISesana et al.ll2007b|) . there are 
some SMBH assembly scenarios that may result in far fewer LISA events (i.e. i f seeds are rare and 



grow primarily through r apid accretion or are very massive already at formation : lTanaka & Haiman 



20091 : iLippai et al.ll2009|) . It is therefore plausible that EM surveys could deliver a larger SMBH 
binary sample than available from GWs. Third, several transient EM surveys are already under 
way, or are planned to be completed before the expected launch date of the LISA mission around 
2020. 

If luminous quasar activity is triggered by major mergers of galaxies, as argued above, then 
the viscous afterglow could plausibly be interpreted as the signature of the birth of a quasar. In this 
paper, we estimate the number of identifiable afte rglow sources, i.e. birthi ng quasars , in the sky, by 
(i) adopting an idealized time-dependent model (iTanaka & Menoull2010l hereafter ITM10|) of the 
evolution of the disk structure, to calculate photometric light curve and variability of the afterglow, 
and (ii) by using the observed luminosity function of quasars as a proxy for the SMBHB merger 
rate. Our two main goals are: 



1 . To assess whether there is any hope of detecting and identifying the viscous afterglows with 
conventional EM telescopes alone. 

2. To see how the identifiability of the afterglows depends on theoretical parameters and to 
delineate the ideal survey attributes (wavelength, angular coverage and depth). We compare 
the derived attributes to those similar to planned large surveys of the transient sky: a soft-X- 
ray survey with specs similar to those that were proposed recently, unsuccessfully, for the the 
Lobster-Eye Wide-Field X-ray Telescope^ (LWFT) mission; and the Large Synoptic Survey 
Telescope^ (LSST) in the optical. 



We find that the detectability of the afterglow is sensitive to the properties of the circumbinary 
disk, in particular to the ratio of the viscous stress to the gas pressure, and to the surface density of 
the disk. We conclude that purely EM identification of the afterglows by the planned surveys are 
unlikely, unless the surface density and the viscosity in the circumbinary disk are at the high end 
of the expected range. In this latter, optimistic scenario, several dozen birthing quasars could be 
identified in a soft X-ray transient survey. We also find that if ^ 1 % of the X-ray radiation emitted 
in the central regions is reprocessed into the optical frequencies by dust surrounding the source. 



^http : //w w w. s tar. le . ac . uk/1 wf t/ 
%ttp://www.lsst.org/lsst 



-5- 



or by warps or geometric splaying in the disk itself (|TM10|) . several dozen afterglows could be 
detected in an optical transient survey, such as LSST. 

This paper is organized as follows. In § 2, we summarize the viscous afterglow model, and de- 
scribe our methods for estimating the identifiable population of AGN harboring a recently merged 
SMBHB. In § 3, we present estimates for the number of identifiable afterglow sources in the sky. 
We summarize our results and offer our conclusions in § 4. 



2. A Simple Model for the Afterglow Population 



In this section, we describe the model and underlying assumptions used to estimate the num- 
ber of observable a fterglow sources in the sky. We use the idealized Newtonian viscous evolution 
model described by iTMld to calculate the light-curve and s pectral evolution of th e source, and to 
obtain light curves in fixed frequency bands. Then, following iHaiman et al.l (|2009h we assume that 
every SMBHB merger ultimately leads to a b right quasar episode. Using the bolometric quasar 
luminosity function of iHopkins et al.l (|2007b|) as a proxy for the underlying SMBH mass func- 
tion, we estimate the number of afterglow sources in the sky that exhibit identifiable photometric 
brightening. 

Throughout this paper, M = Mi +M2, q = M1/M2 and a denote the binary mass, mass ratio and 
semi-major axis, respectively. The physical constants G, c, and have their usual meanings. 



2.1. Modeling Afterglow Light Curves 



Below, we brie fly recapitulate the mai n featu res of th e viscous afterg low model, which was 
advanced by ImPOSI and elaborated upon by ItMIOI (see also lShapirdl2010|) . We refer the reader to 
those earlier works for further details; a derivation of the disk evolutionary equations, in particular, 
can be found in Appendix A of lTMlO . 



A SMBHB in a circular orbit of semi-major axis a will open an annu lar gap in a thin circumbi - 
nary disk at a radial distance i?waii ~ 2a from the binary's center of mass (lArtymowicz et al.lll99l|) . 
The binary will shrink gradually by depositing orbital angular momentum in the disk, maintaining 



^ As discussed bv lHaiman et al.l (l2009h . an alternative method to construct the luminosity function is to calculate 
the SMBH merger rate from the merger rate of dark matter halos, and apply the post-merger light curve L(t,M, q) to 
each merger event. That method requires a third ingredient, namely a way to prescribe the SMBH masses from the 
halo masses. Given the approximate nature of our calculation, we choose the simpler method described in the text. 
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a nearly self-similar geometry with i?waii/(2a) ~ 1. The kinematic viscosity z/gas of the gas com- 
prising the disk is a weak function of radius R, and thus the viscous time ?visc(^) = (2/3)7?^/z^gas 
roughly scales just outside the cavity as ?visc(^waii) oc a^. Because GW emission shrinks the orbit on 
a timescale few = {dlna/dt)~^ oc a"^, there exists a value of a inside which the binary orbit begins 
to shrink faster than the circumbinary gas can viscously follow. This critical binary separation is 
of the order a ~ lOOGM/c^; past this point, the evolution of the binary and the disk are decoupled. 

After decoupling, the evolution of the disk su rface density S can be described by the standard 
equation for a thin, Keplerian accretion disk (e.g., Frank et al.ll2002|) : 



dt ^ ' RdR 



OR ^ ' 



(1) 



In the special case when the dependence of the gas viscosity on radius can be approximated as 
a power-law, v oc R", the su rface density evolution E(i?,0 can b e solved semi-analytically with a 
Green's function formalism (iLynden-Bell & Pringlell 19741 iTMlOl) . starting from an arbitrary initial 
distribution = 0). 

The subsequent evolution of the circumbinary disk is sensitive to several properties of the 
binary-plus-disk system. For a given total binary mass M = Mg x IO^Mq, the relevant observable 
quantities can be expressed in terms the following eight system parameters. 

• The parameter ( = 4^(1-1-^)"^ < 1 is the symmetric mass ratio of the binary, scaled to unity 
for equal-mass binaries. 

• The ratio of the diameters of the circumbinary cavity and the binary separation at decoupling 
is specified by the parameter A = 7?waii/(2a) ~ 1. 

• The gas viscosity is parametrized via the ra tio of the viscous stress to the gas pressure, 
a = i^gas^ksT / Pgas ( Shakura & Sunyaevlll973 ). The choice to scale the viscosity with the gas 
pressure and not the total pressure is motivated by calc ulations that suggest that radiation 



press ure-dominated thin disks may be thermally unstable ( Shakura & Sunyaevlll976 : Pringlc 
19761)3 



The parameter f3 gives the ratio at decoupling of the timescale tow on which the binary 
shrinks due to GW emission to the viscous time ?visc of the gas at the cavity wall. ImPOS 



prescribed /3 = 0.1, based on the fact that the gas near the cavity wall has a very steep den- 
sity gradient and will move on a characteristic time-scale ~ 0. l^visc after the binary torques 



^ Even if thermally stable ( Hirose et al. 20091) . such disks may still be viscously unstable ( Li^htman & Eardlev 

197ilPiraJl9"7ilHirose et Zll2009l) . 
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vanish. However, the calculations by TMld suggest that the lower value of /? ~ 0.05 i 



IS more 

appropriate, as it gives a closer approximation to the time elapsing between the merger and 
when the binary torques cease to influence the gas. 

The disk porosity parameter, 9, prescribes how the optical thickness r between the disk 
midplane and the height where the emitted photons are thermalized relates to the surface 
density: r = 9T.aes, where ergs is the electron- scattering cross section. 

• The viscous evolution model assumes a radial power-law for the kinematic viscosity, z/gas oc 
i?". Prior to merger, viscosity in a circumbinary decretion disk may be expected to be a weak 
function of radius , with n ^0.1; after the merger, n ~ 0.4 may be expected in the accretion 



region of interest (ITM10|) 



Another parameter of interest is the radial power-law index m = d(lnuT,)/dlnR of the disk 
mass profile at decoupling. The standard steady-state solution for a thin accretion disk 
around a single central object satisfies m = 0, with mass flow M = Svrz/gasS constant with ra- 
dius. However, at the cavity wall the binary torques decrete the gas outward; in this regime, 
the disk woul d more likely satisfy m = -1 /2, with torque density z/gasSi?^^^ = constant (e.g.. 



Pringlelll99l|) . We prescribe initial surface density profiles satisfying 

S(^^^wall) < S(i?>^wall) CX RT 



-n+m 



with a steep exponential drop-off in S near R ^ i?waii (|MacFadyen & Milosavljevid 12008 : 



TMIO ). 



• The value of the surface density Ewaii at the cavity wall at decoupling depends on how much 
gas has piled up due to accretion of outer gas and decretion of the gas just inside. We 
parametrize this value as E^aii = 5MHdd/(3vrz/), where Medd is the accretion rate that would 
produce a luminosity corresponding to the Eddington limit for the binary mass M, assuming 
a radiative efficiency of ?] = 10%. The parameter S can be thought of as the product of two 
quantities: the mass supply rate M to the disk in Eddington units, and the enhancement of 
the disk surface density just outside the cavity due to mass accumulation. The second factor 
may be expected to exceed unity. 

In addition to the above parameters, the cosmological redshift z of the source will also obvi- 
ously affect the observed spectra, luminosity, and brightening rates of the sources. 

For a given set of parameters, we use the corresponding surface density evolution S(i?,0 to 
calculate the disk temperature profile at depths where the emitted photons are thermalized, which 
in turn is used to obtain the monochromatic luminosity in the source rest frame: 

L^(0~2x / 71 — —B^l-KRdR. (2) 
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The leading factor of two on the right hand side of equation ^ accounts for the fact that the 
disk radiates from two faces; is the Planck function; €,y = Kabs,i//('«abs,y + «:es) is the ratio of the 
absorption to the total opacity; and the fraction 2e,J (l + e^) is the deviation of the flux from black- 



body (e.g. jRybicki & Lightmanll 19861 ; lBlaesll2004h . The effective temperature in the inner region 



of the disk is higher than a blackbody disk with the same surface density distribution, and thus it 
produces a harder spectrum. The lower limit of integration, i?isco, is the radius of the innermost 
stable circular orbit, for which we adopt the value 3GM/c^, consistent with the assumption that the 
binary remnant has moderate spin, 0.6 < aspin ^ 0.9. The choice for the upper limit of integration 
is somewhat arbitrary, and does not significantly affect Li,{t) as long as it is sufficiently large; at 
larger radii the flux is significantly lower, and evolves on much longer timescales than the inner 
region originally occupied by the cavity. 

The overall disk evolution timescale is roughly the viscous time at the cavity wall, evaluated 
at decoupling: 

^ 1 on NX i\//l-32 ^0.70 A 2.8^-1.36 o0.70r.-0.98/)-0.34 /ox 

^visc ~ 120yr xMg Q A S ^ ■ (3) 

Although the disk continues to brighten and spectrally harden for ~ ?visc after the merger, the 
most dramatic evolution takes place in the first ~ /3?viscj which corresponds to the faster viscous 
spreading of the sharp density edge dXR < i?waii- 

The time-dependent spectrum of the viscously spreading disk can be divided into two parts: a 
nearly static low-frequency component produced predominantly by gas outside i?waii; and a rapidly 
evolving high-frequency component produced predominantly by the gas flowing into the central 
cavity. The characteristic frequency that marks the boundary between the static and variable parts 



of the spectrum can be approximated in the source rest frame (|TM10|) as 

hp,^ ~ 30 eV X Mg-o-^^C-o-^^A^-^a^-f , (4) 
or, in terms of the wavelength. 



Avar ~ 40 nm X M^'X >^ cTi S'""- ■ (5) 

Above, we have defined a-i = a/0.1, = /3/0.1, and 6q2 = 9/0.2. These equations already 
reveal that significant brightening will occur primarily at photon energies in the hard UV to X- 
ray bands. While the characteristic frequency in eq. © can move into the optical bands for the 
most massive SMBHs, the overall evolutionary timescale for these sources, eq. (|3]), becomes long. 
The most rapid evolution takes place as the cavity fills; once it is filled, the system evolves more 
gradually to approach a standard steady thin accretion disk solution around a single SMBH, with 
the quantity z^L,^ peaking at a frequency of roughly t'peak ~ 15//var- The spectrum falls off steeply at 
higher frequencies, and is likely unobservable above u > Sup^ak- 
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Prior to decoupling, the luminosity of the disk at frequencies p > z^var is negligible. Once the 
cavity is filled, the monochromatic luminosity at frequencies z/yar ^ z^ ^ t'peak reaches 

VL, ~ 5-30 X 10^2 e,g,-1^0.92^-042 ^-1.7^0.34^-042^1.2^0.^^^ 

In order to be identifiable in a survey, an afterglow source must exhibit significant brightening, at 
least comparable to the typical variability of typical AGN, at luminosities and frequencies to which 
the survey is sensitive. 



2.2. Modeling the Population of Afterglow Sources 



We now turn to estimating the number of identifiable a f terglo w sources. The approach de- 
scribed below closely follows that described in lHaiman et al.l (|2009l) for estimating the number of 
pre-merger sources that may be detectable by their periodic variability. We begin by prescribing 
the quasar luminosity function as a pr oxy for the SMBH ma ss function. Specifically, we adopt 
the bolometric luminosity function of iHopkins et al.l (|2007bh . and suppose that during a typical 
bright quasar phase, the luminosity and SMBH mass are related via a simple approximate relation, 
L{M) ~ fEAdJ-^-EAdiM), where L-Edd{M) is the Eddington luminosity for an object with mass M, and 
/edd is a constant. This is an admittedly rough estimate, as f^AA is known to have a non-negligible 
spread among the population of bright quasars. However, our simple estimate is sufficient for a 
proof-of-concept; a more precise calculation is not warranted, given the approximate and highly 
idealized nature of the afterglow model a nd the uncertainty in t he system parameters. We choose 
the fiducial values f-^AA ~ 5 ~ 0.3 (e-g-, KoUmeier et al.ll2006|) and assume a rest-frame quasar 
lifetime of ?q ~ 10^ yr (e.g.. lMartinill2004 . 



We further assume that there is a one-to-one correspondence between a SMBHB merger and 
a bright quasar episode, i.e. that a SMBHB merger ultimately triggers quasar activity. This as- 
sumption is consistent with our post-merger disk evolution model, which naturally leads to a state 
with a fully extended disk around a single SMBH, as long as fuel remains available to maintain 
near-Eddington accretion rate at the outer edge of the disk. Given the comoving number density of 
AGN dnp^cn/dM, we are interested in the subset of SMBH merger remnants that are producing an 
observable, brightening afterglow, and have not yet settled to a later, presumably steadier quasar 
phase. To estimate this fraction of SMBHs, we use the afterglow model described above, and cal- 
culate the duration fag over which the photometric emission from a circumbinary disk brightens at 
a rate exceeding some threshold value. This threshold should be chosen to correspond to a bright- 
ening rate that is not only measurable, but is also distinguishable from other possible sources of 
time-variability. The number of variable sources A^ag in the sky of such SMBHB remnants in a 
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given mass and redshift range is then 



dflAGN 

dlnM 



dlnM, 



(7) 



where AV(zio,Zhi) is the cosmological comoving volume between redshifts z\o and Zhi and hagn is 
the space density of SMBHs of mass M, inferred from the quasar luminosity function evaluated at 
L{M). It is worth emphasizing that this expression does not assume that the birthing quasars have 
a bolometric luminosity of L{M) - rather, L(M) here represents the characteristic luminosity that 
is produced by SMBHs of mass M in their quasi-steady quasar state; this asymptotic luminosity is 
reached only > 100 yrs after the merger, according to our afterglow models. 



3. Results and Discussion 

3.1. Basic Parameter Dependencies 

The dependence of the number of detectable variable sources on the various model parame- 
ters for the binary-plus-disk population is of obvious interest, and is not trivial, as each parameter 
affects differently the luminosity, spectral frequency and brightening rate of the afterglow. For ex- 
ample, increasing the binary mass M increases the source luminosity and lowers the characteristic 
frequency of the source, while making the afterglow evolve more slowly — thus, the brightening 
and hardening rates both decrease (making identification more difficult) while the total flux and 
the total afterglow lifetime both increase (making a detection easier). 

To illustrate how each of the parameters and the source redshift affect the detectability of 
variable afterglow sources, in Figure [T] we first plot the basic quantity ?ag.obs = il+z)tag, represent- 
ing the amount of time sources are observed to spend at or above the required threshold for the 
brightening rate. The threshold in this figure is set at JlnLx/(i?obs of at least 10% yr"\ in the soft 
X-ray frequency window of 0. 1 -3.5 keV (motivated by the proposed all-sky X-ray transience sur- 
vey LWFT; see below). The solid black curve in each panel shows ?ag,obs for the fiducial parameter 
values q = a_i = 60,2 = S = X = 1, 13 = 0.05, n = 0.4, m = -1/2, and z = 2. In each panel, we show how 
the apparent duration of the rapidly brightening phase is affected by changes (dashed and dotted 
lines) in one of the system parameters. 

In Figure [21 we plot the corresponding number dN^g/dlnM of sources that exhibit a band 
luminosity Lx of at least 10"^" ergs"^ and an observed brightening rate dlnLx/dtobs of at least 
10% yr"^ in the same 0.1 -3.5 keV frequency window. This is given by the product of the duty 
cycle ?ag,obs/?e space density of SMBHs (eq. |7]), except that a further cut is imposed in 

luminosity. This is because the brightening rate initially may exceed the threshold for SMBHs with 
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masses near the low-mass end of the range shown in the figure, but their band fluxes are still below 
the imposed threshold; these sources are excluded by subtracting the duration of this initial, sub- 
luminous state from the duty cycle ?ag,obs/^G- The number is computed for the whole sky, and over 
a redshift range 1 < z < 3. The line-style scheme is the same as in Figure[T} the solid black curves 
show dNag/dlnM for fiducial parameter values, and the dashed and dotted curves show how the 
number of rapidly brightening sources depends on each parameter. The optimistic luminosity and 
brightening thresholds in this figure are chosen purely for demonstrative purposes. For reference, 
a source with Lx ~ 10^^° ergs"^ and z ~ 1 would, in fact, have a flux of only < 10~'^ ergs~^ cm"^, 
and thus be too faint to be monitored for variability. Also, AGN have been observed to vary in 
their X-ray brightness by as much as order unity on timescales of years. Although the afterglows 
in question would exhibit a monotonic increase in X-ray brightness, along with a corresponding 
monotonic spectral hardening, it is unclear whether a brightening rate of 10% yr"\ even if sustained 
for several years, and accompanied by a monotonic hardening of the spectrum, would be sufficient 
to distinguish an afterglow candidate from other X-ray variable sources. 

The steep cutoff at high binary masses seen in both tag and dNag/dlnM, in Figures [T] and [21 
respectively, has two causes. One is that for sufficiently large M, the emission frequency of the 
source becomes too low, and falls out of the soft X-ray window. The other reason is that the disk 
evolution timescale t^i^^ scales as oc M^^, so that for sufficiently large M the disk evolves so slowly 
that its brightening rate never reaches the adopted threshold of 10% yr"' . 

Figures [T] and [2] also show that the duration of the brightening phase, and the mass function 
of the afterglow sources depend strongly only on the parameters a, (3 and S. This is due to the 
fact that the afterglow frequency range and evolution timescale scale steeply with these parameters 
(see equations [3] and H]). Increasing a, increasing S and decreasing (3 relative to their fiducial 
values all have the effect of increasing the afterglow emission frequencies and pushing it further 
into the survey frequency window, while also increasing the brightening rate of the afterglow. 
The quantities z/var and ?visc both depend only weakly on 9, and ( varies too weakly in the range 
0. 1 ^ ^ < 1 to have a sizable effect. Increasing (decreasing) the parameter A results in afterglows 
that are further inside (outside) the frequency window but evolve much more slowly (quickly) - 
the two effects tend to cancel out, and yield a relatively weak overall effect on dNag/dlnM. 

It is worth noting that the brightest afterglow sources satisfying a fixed Jin Lx/<i?obs > 10% yr"^ 
are not the most massive ones. This is because disks around more massive BHs evolve more slowly, 
and the brightening rate is greatest early in the post-decoupling disk evolution when the source is 
dimmer. For most of the parameter value combinations probed in Figure[2] — excepting a, f3 and S 
for the moment — the most luminous sources brightening at or above the threshold rate are those 
with binary masses of (0.5-2) x IO^Mq. Interestingly, this mass range lies in the middle of LISA'S 
sensitivity window. 
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The maximum band luminosities of these sources are typically 1 - 4 x 10"^^ erg s"^ , and behave 
roughly as described in equation Q. For a = I and 5 = 3, the masses and luminosities for the 
brightest afterglow sources are somewhat greater: > 5 x 10^ Mq and > 10"^"^ ergs"^ 

The range of luminosities across the parameter combinations probed in Figures [Hand [21 corre- 
spond to an observed soft X-ray flux of Fx ~ 10"^^- 10"'"^ x S^-^ erg s"^ cm"^ in the range 1 ^ z ;^ 3. 
Thus, if the approximate location of the source is known via a prior GW detection, the afterglow 
would be observable during the rapidly brightening phase at the sensitivity achieved by existing 
instruments such as XMM -Newton, ROSAT HRI, and Chandra (at ~ 100 ks exposure; see, e.g.. 



Brandt & Hasingenl2005h . For the parameter combinations probed in Figures [T] and [2l there are 



at most ~ 100 sources in the sky with Lx > 10"*^ ergs ' and dlnLx/dtobs > 10% yr \ and they 
maintain this luminosity and brightening rate for ?ag,obs ^ 10 yr in the observer frame. 



3.2. Counts of Birthing Quasars in X-ray and Optical Surveys 

Applying the methods described above to estimate the afterglow light-curve and source popu- 
lation, we next calculate the source counts of identifiable afterglows as a function of their apparent 
flux. Based on the fact that the mergers of massive dark matter halos peak at a redshift z ~ 2, 
we limit our analysis to AGN in the redshift range 1 < z < 3. We assume that the typical disk- 
plus-binary system has parameter values ^ = 0. 1, (3 = 0.05, X= I, 6 = 0.2 n = 0.4, and m = -1/2. 
In reality, these parameters will vary from system to system, perhaps by a great deal. However, 
because our calculations above (equation [6] and Figure [2l) suggest that neither the luminosity nor 
the mass function of afterglow sources are likely to be affected by the values of these parameters 
by more than an order of magnitude, we limit further exploration of the parameter space to the a-S 
plane. The value of the viscosity parameter a is highly uncertain; however, numerical simulations 
of MHD disks suggest in the radiation-dominated r egions of an accre tion disk a ~ 1 is consis- 



tent with our type of viscosity prescription u oc Pg^s ([Hirose et al.ll2009|) . Assuming that the mass 
supply rate around the afterglow phase is at least comparable to the rate during prolonged AGN 
activity, we expect S to be at least as great as the ratio of the typical luminosity of the typical AGN 
episode to the Eddington luminosity limit of the SMBH engine, i.e. S > /ndd = 0.3. Depending on 
how much mass the circumbinary disk accumulates near its inner wall prior to decoupling, S may, 
however, significantly exceed /edd- 
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3.2.1. X-ray Surveys 

We consider a hypothetical survey with observational capabilities similar to those of the 
proposed LWFT mission, sensitive to photon energies of 0.1-3.5 keV, down to fluxes of ~ 
3.5x10"'"* ergs"' cm"^ (confusion limit for LWFT). 

In Figure |3l we present the X-ray source counts of objects increasing in band luminosity by 
at least 3% yr"\ 10% yr"', 30% yr"' and 100% yr"' in the relevant frequency band. (Note that 
the energy range approximately coincides with the sensitivity window of existing deeper-exposure 
telescopes such as Chandra and XMM-Newton, and the planned all-sky survey eROSITA.) In the 
figure, the height of the histogram pillars are the number of sources in each logarithmic flux bin, 
with each bin having a width AlogjQ(F erg"' cm^s) = 0.5. The dashed vertical line in each panel 
shows the approximate source confusion flux limit. These results suggest that only if S is large, i.e. 
if there is significant mass accumulation at the inner wall of the circumbinary disk, then rapidly 
brightening afterglow sources could be detected by the model mission: several sources at 5 ~ 3 
and as many as dozens of sources for S ~ 10. The brightest of these sources have central SMBH 
masses of M > 1O''5' ^M0, and thus most of the X -ray-detectable sources would also be observable 
by LISA. However, for 5* > 3, some would fall outside USA's sensitivity window. These sources are 
expected to continue to brighten at a slightly reduced rate at harder frequencies, hv > 3.5 keV; this 
is a prediction that could be tested with pointed follow-up observations. Our calculations indicate 
that birthing quasars will be difficult to identify with existing and planned wide-angle soft-X-ray 
surveys. For example, the eROSITA all-sky survey is expected to have a semi-annual flux limit of 
> 10"'^ ergs"'cm"^ with a time resolution of ~ months at those flux levels. Intrinsic absorption of 
soft X-rays by the bir thing quasar itself could also be an observational barrier for at least some of 



the sources (see, e.g. jBrandt et al.ll2000l) . especially if galaxies harboring merging SMBHBs tend 



to be more heavily shrouded in gas and dust than the general population of active galaxies. 

Whether the monotonic brightening of the afterglow would be sufficiently distinguishable 
from the X-ray variability of faint AGN and other sources is an open (and more difficult) ques- 
tion. Luminous X-ray AGN have been known to vary by as much as ~ 10- 100% yr"' (e.g., 
Mushotzky et al.|[l993 l). Any monotonic brightening in the X-ray must be distinguished from other 



sources of intrinsic variability, in addition to any instrumental error close to the faint-flux detection 
limit. It would appear likely, however, that a source that monotonically increases its X-ray lumi- 
nosity by up to order unity per year for several years, while showing a consistent and monotonic 
hardening of its spectrum for the entire duration, would be fairly unusual. While an increase in 
the accretion rate through a standard thin accretion disk is also expected to produce a simultaneous 
brightening and spectral hardening, the spectral evolution will be different in detail. We also note 
that a variability survey , utilizing the ROSAT all sky survey, as well as pointed ROSAT observa- 



tions (|Grupe et al.ll200l|) found no correlation between changes in flux and spectral hardening in a 
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sample of luminous soft-X-ray AGN. If the brightening is caused by the viscous afterglow, it would 
also slow down on humanly tractable timescales and the light curve could be checked at different 
observational frequencies against the evolutionary models for the viscously spreading disk. This 
feature would also help distinguish birthing quasars from tidal disruption events. 

It is worth emphasizing that a major caveat of the above analysis is that the thin-disk formalism 
adopted in our afterglow model breaks down for models with 5 > 3. Indeed, in this massive 
regime, the disk midplane temperature becomes sufficiently high that a one-dimensional estimate 
for the disk scale-height-to-radius ratio H /R, evaluated at the cavity wall, exceeds unity (iMPOSi : 
TMlOh . This suggests that horizontal advection would become a significant factor in determining 
the structure of the inner accretion flow. The evolution and emission properties of the disk in this 
advective regime is highly uncertain and it remains a subject of active resea rch. Advective disks 
have generally been associated with radiatively inefficient accretion states (see lNarayan & Ouataert 
2005L for a review), but the situation may be different h ere, wit h a high density accretion flow and 
an intrinsically time-dependent flow. A s suggested by ImPOSI. ho rizontal advection could act to 
make the disk "slim", as in the models of Abramowicz et al. and thus relatively radiatively 



efficient. It is also unclear how well binary torques can open a central cavity in the first place, if the 
disk becomes geometrically thick. Adding to the various modeling uncertainties is the possibility 
that a sizable fraction of the viscously dissipated energy could be deposited into a hot disk corona 
(e.g.. iLiu et al.ll2003a|) . which would lead to a significant increase of the afterglow high energy 
emission, well above what is predicted by our strictly thermal afterglow models. 



3.2.2. Optical Surveys (LSST) 

Our simple viscous emission model predicts that the same SMBH binary remnant that pro- 
duces an X-ray afterglow would brighten at lowe r, optical frequencies at a rate of > a few% yr"^ 



several years prior to the X-ray afterglow ( TMlot) . This less pronounced varia bility is comparabl e 



in magnitude to the intrinsic r.m.s. long-term variability of optical AGN (e.g. JUlrich et al.lll997h . 
However, if the approximate location of the source is known through a GW detection, searching 
for AGN exhibiting steady optical/ultraviolet brightening may identify the source before the X-ray 
afterglow. Alternatively, in the absence of a GW signal, a wide-angle variability survey of optical 
AGN with high time resolution, such as those possible with LSST, could possibly still be used to 
select afterglow candidates for follow-up X-ray observations. 

Motivated by these possibilities, we apply our simple source-count estimate to the LSST u 
photometric band (330-400 nm). We choose this band because it is the LSST filter with the 
highest frequency range, and thus the one in which the afterglow brightening is likely to be the 
most prominent. The results are presented in Figure IH In both panels, the dashed vertical lines 
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demarcate the limiting flux for signal-to-noise of 50 over a ~ 1 yr period (~ 450s accumulated 
exposure). With the fiducial parameter choices in the left panel, there is clearly no hope of a 
detection. We also find that, unlike in the X-ray bands, raising the value of S does not increase the 
source counts. This is because while increasing S beyond the fiducial value pushes the evolving 
portion of the spectrum into the model instrument band, it pushes it further out of the LSST u 
band, reducing the optical variability. Therefore, we instead show, in the right panel, the expected 
number counts for the larger value « = 1, on which the emission frequency depends less strongly. 
We find that the sources are still likely to be too dim and too few to be identifiable with high 
confidence from among the large number of AGN expected to be detected by LSST. However, 
gradually brightening optical AGN could still be cataloged, and their X-ray luminosities could 
be cross-checked with data from instruments such as eROSITA for subsequent X-ray afterglows, 
which may still prove useful. 

More promising for optical surveys is the possibility that the X-ray afterglow ma y be pro mptly 



reprocessed into optical or infrared frequencies, either by dust surrounding the source (|MP05|) or by 



warps and vertical splaying in the circumbinary disk (|TM10|) . If a significant portion of the X-ray 
emission is reprocessed, then the afterglow source will appear as an AGN whose optical/infrared 
luminosities can brighten by d\nLjj/dt > 10% yr"^ Such a monotonic variability would exceed 
the typical long-term r.m.s. variability in AGN brightness at these frequencies, and is likely to be 
identifiable by LSST. For purely demonstrative purposes, we adopt here a simple toy model, in 
which a fraction of the emitted power above u > I keV in the source's rest frame is reprocessed 
(thermalized) and re-emitted at frequencies below the ultraviolet energy 10 eV, so that the power 
below hu < 10 eV is enhanced as 

^hu<O.Ql keV ~ ■^'/!i'<0.01 keW + frpLhu>l keV- (8) 

We further assume that the fractional energy enhancement is roughly uniform in the optical and 
infrared — i.e., we ignore line emission from recombination processes — so that for our purposes 
the optical spectral emission is given by 

~ X f 1 -^^->'^ev \ 

V J^hu<O.Ol keV / 



This prescription is similar in spirit to the model of iGierlihski et al.l (|2009|) . who showed that 
reprocessing signature of the X-ray outbursts of the stellar-mass black hole system XTE J1817- 
330 is consistent with reprocessing a constant fraction of the bolometric X-ray luminosity. We 
also neglect the reprocessing time, which is dominated by the light-travel time and much shorter 



than the variability timescales of interest here (see, e.g. JPeterson & Hornell2004|) 



The source counts in the LSST u band for this model, for a reprocessing fraction f^ = \0^, 
are shown in Figure [H Our simple calculations suggest that perhaps dozens of afterglows could be 
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detected if the X-ray emission is reprocessed, for moderately optimistic parameter values, e.g., for 
5 > 1 . It is worth cautioning that the reprocessed fraction f^p is highly dependent on th e vertic al disk 
geometry, which itself may be rapidly evolving during the afterglow. For instance, ItMIO found 
in their non-irradiated afterglow models that the scale-height-to-radius ratio H/R can be a steeply 
increasing function of radius during the period when the cavity is refilling, but not necessarily 
before or after this phase. Such complications should be included, along with details of radiative 
transfer, in a more realistic analysis of disk irradiation and reprocessing. Finally, absorption of the 
reprocessed UV/optical emission by gas and dust, surrounding the nuclear SMBH on larger scales, 
could be another observational hindrance for at least some of the sources. 



4. Conclusions 

Using an idealized model for the population of coalescing SMBHBs, and for the light curve 
of the afterglow produced by the viscously spreading post-merger circumbinary disk, we have 
shown that ongoing afterglows of SMBHB mergers may be present in the data sets of wide X-ray 
and optical surveys. In soft X-ray bands, this requires that the surface density and the viscosity 
in the circumbinary disk be at the high end of the expected range, while afterglows could only 
be found in optical surveys if the X-ray emission is promptly and significantly reprocessed into 
optical frequencies. 

Despite the highly approximate nature of our analysis and other model uncertainties, our 
calculations provide a proof-of-concept for a very general hypothesis: SMBHB mergers may exhibit 
identifiable, steady brightening rate for a period of the order of decades, and such afterglows could 
be detected serendipitously in a large survey that revisits the sky at least every few months for 
several years. Our more specific findings can be summarized as follows: 

• For optimistic parameter values, several birthing quasars, brightening by at least d InLx/dtob^ > 
30% yr"' for several years, could be identified in the 0. 1 - 3.5 keV soft X-ray band by an all- 
sky survey with specifications comparable to those proposed for the LWFT mission. 

• At any given time, there could be up to Ni,g ~ 100 sources in the sky that exhibit a brightening 
at or above dlnLx/dtohs > 10% yr"' , with soft X-ray luminosities Lx > 10"^^ erg s"^ . The most 
luminous sources typically spend ?ag.obs ^ 10 yr in this state, and thus can be monitored on 
humanly tractable timescales. These numbers depend weakly on most system parameters. 

• To have any hope of detecting birthing quasars, a survey has to reach a depth of at least a few 
X 10"^^ ergs"' cm"^. However, the slopes of our calculated logA'^-logS distributions at fluxes 
just below this threshold are relatively shallow (Figures [31 H] and [5]), implying that surveys 
should favor large angular sky coverage over depth, once they reach this flux threshold. 
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• If identified, candidate sources can be followed up by pointed observations at higher fre- 
quencies, where they are expected to continue both their monotonic brightening and their 
spectral hardening. 

• Most birthing quasars that are identifiable have, coincidentally, SMBH masses lying in the 
middle of USA's sensitivity window (M ~ IO^Mq), and are thus members of the same popu- 
lation that would be probed with GW detections. However, a minority (> few % for 5 > 3) of 
the detectable X-ray variables have masses of ^ 10^ Mq, probing a population above USA's 
range. 



• These sources may be identifiable by LSST if a fraction as low as 
promptly reprocessed into the optical frequencies. 



1% of the X-ray flux is 



Our calculations are contingent on theoretical caveats of the afterglow scenario we have con- 
sidered. The two primary uncertainties regarding the post-merger evolution of the circumbinary 
cavity are related to the viscous and advective properties of the disk. As stated in § 12.11 the viscos- 
ity of accretion flows, including the possibility of viscous instability, are not well understood when 
radiation pressure dominates gas pressure, which is the relevant regim e for the gas refi lling the cir- 
cumbinary cavity. Additionally, the disk may be geometrically thick (|MP05l : ItM10|) . either right 
at decoupling or later during the afterglow phase, suggesting that horizontal advection may play a 
significant role in determining the surface density evolution and the disk net emission properties. 
The importance of viscous instabilities in radiation-dominated accretion flows remains a general 
open question, and the role of advection in a viscously spreading accretion flow remains a largely 
unexplored regime. More detailed studies of the circumbinary cavity will be needed to address how 
these effects may affect the emission predicted by simple analyses based on a thin disk formalism 
such as ours. Another major uncertainty is the validity of our assumption that quasar activity can 
be associated with SMBH coalescence. In reality, there may not be a one-to-one relation: it is 
possible that for at least some AGN, gas accretion or changes in radiative efficiency are triggered 
by mechanisms other than SMBH mergers; conversely, some SMBH mergers may not trigger pro- 
longed quasar activity. If the former is true, our analysis overestimates the number of identifiable 
afterglow sources; if the latter is true, then our results could in principle be an underestimate. 

For completeness, we note that while we focused here on the viscous afterglows, other SMBHB 
merger-related signatures could also be looked for in EM surveys. For example, the GW-emission- 
induced mass-loss and recoil c an cause strong disturbances in the circumbinary disk, which can 



produce a detectable afterglow (Lippai et al. 20081; Ischni ttman & Krolik 2008 ; Shields & Bonning 



2009; Corrales et al. 



201 



I 



Rossi et al 



20101) . For the 



20081 : lO'NeiU et al.ll2009l : iMegevand et al 
low SMBH masses of ~ lO ^MfD relevant for LIS A, these signatures are expected to have a short 
duration ~ few years (e.g., ICorrales et al.l 120101) and would be too rare to be found serendipi - 



tously, without a trigger from LISA. However lSchnittman & KrolikI (|2008|) and lShields & Bonning 
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(|2008|) focused on these signatures in disks around more massive SMBHs, which occur on longer 
(~ lO'^yr) time-scales, and proposed detecting a flare by monitoring a population of AGN in the 
infrared or X-rays bands. Another possibility is that the binar y is activated, a nd pro duces periodic 
emission, tracking the orbital frequency, prior to the merger. iHaiman et al. argued that as 



long as this emission is at a few percent of the Eddington luminosity, a population of these variable 
sources, with periods of tens of weeks, may be identifiable in optical or X-ray surveys. 

To conclude, the concomitant observation of a SMBHB merger based on GW and EM signals 
remains by far the most promising scenario for the unambiguous detection of such systems. The 
precision with which LISA would determine the masses, spins, and luminosity distances of coa- 
lescing binaries can not be replicated by current or planned EM telescopes. However, detections 
based on EM signatures alone could still help identify SMBHB mergers before LISA is launched, 
and perhaps more importantly, possibly outside LISA'S mass sensitivity window. Detecting the 
EM signatures from the mergers of the most massive SMBHs would complement the synergistic 
EM-plus-GW observations of lower-mass systems, and help provide a more complete picture of 
the accretion physics and cosmological evolution history of SMBHBs. 
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Fig. 1. — The approximate amount of time spent by an afterglow source with an observable band 
brightening rate of at least 10% yr"^ in the 0. 1-3.5 keV photon energy range. The most luminous 
such sources spend t^g > \0 yr in this brightening phase, roughly independently of the model 
parameter values. The solid curve in each panel shows the fiducial set of parameters, q = a_i = A = 
S = 6*0.2 = !,/? = 0.05, z = 2, n = 0.4 and m = -1/2. The dotted and dashed curves show how the 
duration of the rapidly brightening phase is affected by changes in each of the model parameters 
(see § 2.1 for a definition of each parameter). 
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Fig. 2. — The approximate number of afterglow sources in the sky at any given moment whose 
luminosity in the 0.1-3.5 keV energy range is (1) at l east 10"^° ergs" ^ and (2 ) increasing by at 
least 10% yr"' . We use the AGN luminosity function of Hopkins et al. (2007W) as a proxy for the 
SMBH mass function, and associate each episodic activity of AGN with a SMBH me rger whose 
afterglow light curve is given by the time-dependent model of iTanaka & Menoul (|2010f) . The solid 
curve in each panel shows the fiducial set of parameters, q = a-i = X = S = 6*0.2 = ^, P = 0.05, 
1.5 < z < 2.5, n = OA and m = -\/2. The dotted and dashed curves show how the number of sources 
are affected by changes in each of the parameters. Note that the most massive sources are not the 
most luminous in the band of interest (see text). For plausible parameter values, the most luminous 
sources have masses, band luminosities, and number in the sky of > IO^Mq, L > 10"*^ ergs"^ and 
overall number A'^ag ^ 100, respectively. The masses of luminous sources of interest coincide with 
the LISA sensitivity window, 10^ - 10^(1 +z)~^Mq. 
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Fig. 3. — The number of afterglow sources as a function of their 0. 1 - 3.5 keV flux, in the redshift 
range 1 < z < 3. The histogram shows the number of sources in each logarithmic flux bin of width 
AlogjQ(F erg"^ cm^s) = 0.5. We consider sources with IO^Mq < M < IO'^Mq. The histograms 
demarcate, in order of decreasing line thickness, the number of sources in each flux bin exhibiting 
a brightening rate of at least 3% yr"^ 10% yr"', 30% yr~\ and 100% yr"^ The dashed vertical 
line in each panel is the detection limit for point sources in our model survey. All panels have 
parameter values ^ = 0.1, a-i = X = 9o.2 = 1, /5 = 0.05, 1.5 < z < 2.5, n = 0.4 and m = -1/2. We vary 
the value of S, a measure of the disk mass, in each panel: 5 = 1 in panel (a), 5 = 3 in panel (b), and 
5= 10 in panel (c). 
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Fig. 4. — Same as Figure [3] but in the u photometric band (330-400 nm). The dashed vertical 
line is the limiting flux to achieve a signal-to-noise of at least 50 over ~ 1 yr of operation for LSST 
(~ 30 exposures at 15s each). Both panels have parameter values q = 0.\, S = X = ^0.2 = 1, /? = 0.05, 
I < z <3, n = 0.4 and m = -1/2. Panel (a) has a = 0. 1, and panel (b) has a=l (more viscous disk). 
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Fig. 5. — Same as Figure |4] but assuming that /rp = 1% of the power emitted above 1 keV is 
reprocessed into optical and infrared bands (see text for details). The dashed vertical line is the 
limiting flux as in FigurelH Both panels have parameter values q = 0.l, a_i = A = 60,2 = 1, /3 = 0.05, 
I < z < 3, n = 0.4 and m = -1/2. Panel (a) has 5=1, and panel (b) has 5 = 3 (more massive disk). 



